Universality of the Lyapunov regime for the Loschmidt echo 
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The Loschmidt echo (LE) is a magnitude that measures the sensitivity of quantum dynamics to 
perturbations in the Hamiltonian. For a certain regime of the parameters, the LE decays exponen- 
tially with a rate given by the Lyapunov exponent of the underlying classically chaotic system. We 
develop a semiclassical theory, supported by numerical results in a Lorentz gas model, which allows 
us to establish and characterize the universality of this Lyapunov regime. In particular, the univer- 
sality is evidenced by the semiclassical limit of the Fermi wavelength going to zero, the behavior for 
times longer than Ehrenfest time, the insensitivity with respect to the form of the perturbation and 
the behavior of individual (non-averaged) initial conditions. Finally, by elaborating a semiclassical 
approximation to the Wigner function, we are able to distinguish between classical and quantum 
origin for the different terms of the LE. This approach renders an understanding for the persistence 
of the Lyapunov regime after the Ehrenfest time, as well as a reinterpretation of our results in terms 
of the quantum-classical transition. 

PACS numbers: PACS: 03.65.Sq; 05.45.+b; 05.45.Mt; 03.67.-a 



I. INTRODUCTION 



Controlling the phase in the evolution of a quantum system is a fundamental problem that is becoming increasingly 
relevant in many areas of physics. In relatively simple systems, like a quantum dot in an Aharanov-Bohm ring 1 , the 
phase can even be measured by transport experiments. The development of the quantum information field requires the 
control of the phase of increasingly complex systems 2 . Such a control is hindered by interactions with the environment 
in a way which is not completely understood at present. 

Nuclear Magnetic Resonance provides a privileged framework to test our ideas on the evolution and degradation of 
the quantum phase. The phenomenon of spin echo, through the reversal of the time evolution, allows to study how an 
individual spin, in an ensemble, loses its phase memory 3 . The randomization of its phase appears as a consequence 
of the interaction with other spins that act as an environment. Recently, it has become possible to test the phase 
of the collective many-spin state through the experiments of Magic 4 and Polarization 5 echoes. In these cases a local 
polarization "diffuses" away as consequence of the spin-spin interactions in the effective Hamiltonian H. The whole 
many-body dynamics is then reversed by the sudden transformation TL — » —TL . However, there is an increasing failure 
to reach the initial polarization state which is a consequence of the fluctuations of the phase of the complex quantum 
state 6 and is a measure of the entropy growth 7 . 

Surprisingly, the rate of loss of information of the phase appears as an intrinsic property of the system, being quite 
insensitive to how small is the coupling to the external degrees of freedom or the precision of the reversal 8 . This may 
be interpreted as analogous to the residual resistivity of impure metals. When the direct coupling to the thermal 
bath is decreased by lowering the temperature, the resistivity becomes controlled by the reversible elastic scattering 
with impurities 9,10 . The common feature of both intrinsic behaviors is the complexity of the dynamics that justifies 
the stosszahlansatz or molecular chaos hypothesis. 

However, such an hypothesis does not seem to be compatible with our basic knowledge of quantum dynamics. 
Unlike classical mechanics, quantum dynamics exhibits a remarkable insensitivity to initial conditions 11,12 . That 
is why the field known as quantum chaos deals mainly with the quantum stationary properties of systems whose 
underlying classical dynamics is chaotic. Among these properties, the ones more frequently studied are the level 
statistics 13 , wave function scarring 14 , and parametric correlations 15 . A notable exception among these studies was 
that of Peres 16 , who realized that classically chaotic and integrable systems behave differently under imperfect time 
reversal, for very short and long times. It is through the experimental findings above cited, that the study of time 
evolution of classically chaotic systems has gained a privileged place in nowadays research. 

A simplified version of the echoes experimentally studied is the so-called Loschmidt echo (LE) 
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where | Vo) is an arbitrary initial state that evolves forward in time under the system Hamiltonian Ho for a time t, 
and then backwards under a slightly perturbed Hamiltonian Ho + S. The amplitude m(t) of the LE is the overlap 
between the two slightly different evolutions of the same initial state, and M(t) quantifies the departure from the 
perfect overlap. Because of this important property, within the field of Quantum Information the LE is referred to as 
"fidelity" 17 . Alternatively, M(t) can also be written as the trace of the product of two pure-state density matrices p 
or Wigner functions W evolving with different Hamiltonians, 



M(t) = tr{/^ 0+E (t) p Ho (t)} = {2nh) d J dr J dp W Ho+E (r, p; t) W no (r,p;t) 
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We have used the standard definitions 
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where d is the dimensionality of the space. 

In consistency with the experimental behavior of the polarization echo, the LE of a classically chaotic one-body 
Hamiltonian Ho was found to exhibit an intrinsic decay rate 18 . This result is valid beyond some critical value of the 
perturbation. Interestingly, the decay rate is precisely the Lyapunov exponent A of the classical system. A related 
relevance of the classical dynamics had been hinted from the analysis of the entropy growth of dissipative systems 19 . 

The purely Hamiltonian character of the model of Ref. 18, as well as the result of a classical parameter (A) governing 
a bona fide quantum property (M), attracted considerable attention. A quite intense activity has been devoted in 
the last two years in order to test these predictions in various model systems and pursue further developments of the 
theory 20 " 33 . 

The Lyapunov behavior has been numerically obtained in models of a Lorentz gas 21 , kicked tops 22,23 , Bunimovich 
stadium 24 , bath tube stadium 25 and sawtooth map 26 . The analytical results have been mainly focused in the small 
perturbation region. Jacquod and collaborators 22 identified the regime below the critical perturbation as following 
a Fermi Golden Rule through the energy uncertainty produced by the perturbation which were also analyzed with 
semiclassical tools 27 . Prosen and collaborators 28,29 showed that M(t) in the perturbative regime depends on the 
specific time dependence of the perturbation correlation functions. 

The Lyapunov regime bears a clear signature of the underlying classical dynamics. This observation lead Benenti 
and Casati 26 to propose that the independence of the decay rate on the perturbation strength is a consequence of 
the quantum-classical correspondence principle. As we will analyze in this work, the situation is far less trivial. We 
are going to show that this regime persists for times much larger than the Ehrenfest time (as defined by Berman and 
Zavslavsky 34 ). In addition, the quantum LE is functionally different than what a direct estimation would yield for 
the classical LE (for the chaotic 18 , as well as the integrable 30 cases). Moreover, the classical counterpart of the LE is 
problematic since a wide range of dynamic behaviors is obtained in different situations 31 ' 32 . 

The LE has also been studied in different disordered systems 20,33 . It has been shown in both cases that the long 
range of the perturbing potential, as emphasized in Ref. 18, is crucial in order to obtain a perturbation independent 
regime. 

The various approximations that the semiclassical theory of Ref. 18 relies on were further corroborated using 
an initial momentum representation of the wave-packet 35 . This changes the sum over an uncontrolled number of 
trajectories into only one, which allows the exact numerical evaluation of the semiclassical expression for the echo. 

Taking the perturbation as the action of an external environment allows us to think of the LE as a measure of 
the decoherence. This approach has been advocated by Zurek 36 , and extended 37 by studying the decay of M(t) as 
expressed by a product of Wigner functions (Eq. (1.2)). A semiclassical approximation to the Wigner function allows 
us to separate the different contributions to the LE coming from classical and non-classical processes. As we discuss 
in detail in the sequel, such distinction enables to quantify how decoherence builds in until the classical terms finally 
dominate the LE. 
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With the goal of addressing experimentally relevant systems 3 °, we illustrate our findings in a simple model with 
classical chaotic dynamics: the Lorentz gas. This system has been shown to exhibit a well defined Lyapunov regime 21 . 
The semiclassical theory that we develop, as well as the extensive numerical results that we present in this work, 
allows us to establish and characterize the universality of the perturbation independent regime. 

This universality manifests itself by the robustness of the Lyapunov regime with respect to various effects. Firstly, in 
the semiclassical limit of Fermi wavelength Xf going to zero, the borders of the regime extend from zero perturbation 
up to a classical upper bound. Secondly, and as stated above, for finite Xf the Lyapunov regime extends up to times 
arbitrarily larger than Ehrenfest's time. Finally, universality is also evidenced by the insensitivity of the Lyapunov 
regime with respect to the form of the perturbation or the (non-averaged) behavior of individual semiclassical initial 
conditions. 

The paper is organized as follows: in section II we develop the semiclassical approach to the LE with a quenched 
disorder playing the role of the perturbation, as proposed in Ref. 18. We then discuss the main assumptions and set 
the theoretical framework that will be further developed in the rest of the paper. In section III we consider a specific 
model, the Lorentz gas, and a different perturbation than in the previous case. We first characterize the classical 
dynamics of the Lorentz gas, as well as that of the perturbation, and then present a semiclassical calculation of the 
LE, discussing the different regimes predicted by the theory. In section IV we concentrate the main results of this 
work. The universality of the Lyapunov regime is discussed and supported with numerical results on the semiclassical 
limit, the behavior after the Ehrenfest time and the effects of averaging. In section V we discuss the relation of the 
LE to decoherence by studying the semiclassical approximation to the Wigner function and reinterpreting the results 
of Sec. II under this new highlight. We conclude in section VI with some final remarks. 



II. THE LOSCHMIDT ECHO - SEMICLASSICAL ANALYSIS 



A. Semiclassical evolution 



In this section we calculate the Loschmidt echo (Eq. (1.1)) for a generic chaotic system Ho and a perturbation T, 
arising from a quenched disorder. We follow the analytical scheme of Ref. 18, discussing the main assumptions and 
the generality of the results. We choose as initial state a Gaussian wave- packet (of width a), which is the closest we 
can get to a classical state. 
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We will keep the spatial dimension d arbitrary in the analytical calculations, but it will be fixed to d = 2 for the 
numerical studies of the Sec. IV. It has been shown 41 that if the initial state is a superposition of N Gaussians, 
the final result is the same exponential decay one obtains with a single Gaussian but normalized by N. Thus, the 
assumption of Eq. (2.1) is as general as the decomposition of a given initial state into a sum of Gaussians. The time 
evolution of state ij)(r, 0) is given by 



V(r,i) - J drK(r,r;t) ^(r,0) 



(2.2) 



with the propagator 

K{Y,r;t) = {r\e- mt > h \v). (2.3) 

We will use the semiclassical expansion of the propagator 42,43 as a sum over classical trajectories s(f, r, t), going 
from r to r in a time t, 

K(r,r;t) = ^ K s (r,r;t) , 
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valid in the limit of large energies for which the de Broglie wavelength (Xf — 2ir/kF = 2nh/p ) is the minimal length 
scale. S s (r,r;t) — J Q dtL s (q s (t),q s (t);t) is the action over the trajectory s, and L the Lagrangian. The Jacobian 
C s = |detB s | accounts for the conservation of classical probabilities, with the matrix 
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obtained from the derivatives of the action respect to the various components of the initial and final positions. We 
note fi s the Maslov index, counting the number of conjugate points of the trajectory s. Since we will work with fairly 
concentrated initial wave-packets, we use that Vr i <S , s |r =ro = Ps i (Fs i i s the i-th component of the initial momentum 
of the trajectory s) and we expand the action as 



S s (r, r; t) ~ S s (r, r ; t) - p s ■ (r - r ) 



(2.6) 



We are lead to work with trajectories s that join r to r in a time t, which are slightly modified with respect to the 
original trajectories s(r, r,i). We can therefore write 



^(r,i)= K s(r,r ;t) J dr exp 



s(r ,r,t) 

2\d/i 



-p s • (r-r ) 



^(r,0) 



(4na 2 ) K s (r,r ;t) exp 

s(r ,r,t) 



a I— \ 2 

-^( Ps -Po) 



(2.7) 



where we have neglected second order terms of S in (r — ro ) since we assume that the initial wave packet is much 
larger than the Fermi wavelength (a ^> Xf)- Eq. (2.7) shows that only trajectories with initial momentum p s closer 
than h/a to p are relevant for the propagation of the wave-packet. 

B. Semiclassical Loschmidt echo 

The amplitude of the Loschmidt echo, defined in Eq. (1.1), for the initial condition (2.1), can be approximated 
semiclassically as 
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Without perturbation (S = 0) and restricting ourselves to the terms with s — s (which leaves aside terms with a 
highly oscillating phase) we simply have 
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We have performed the change from the final position variable r to the initial momentum p s using the Jacobian C, 
and then carried out a simple Gaussian integration over the variable p s . 

For perturbations E that are classically weak (as not to change appreciably the trajectories governed by the dynamics 
of Ho), we can also neglect the terms of (2.8) with s^s and write 
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Where AS's is the modification of the action, associated with the trajectory s, by the effect of the perturbation E. It 
can be obtained as 



AS S = - f <ffiS a (q(t),q(t)) , 
Jo 



(2-11) 
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in the case where the perturbation appears as a potential energy in the Hamiltonian (like we discuss in this chapter) . 
If the perturbation is in the kinetic term of the Hamiltonian (like in Sec. Ill), there is an irrelevant change of sign. 

Clearly individual classical trajectories will be exponentially sensitive to perturbations and the diagonal approx- 
imation of Eq. (2.10) would sustain only for logarithmically short times. However, it has been argued 35 that this 
approximation is valid for much longer times because of the structural stability of the manifold 27 which allows for the 
existence of trajectories arriving at r and departing exponentially close to ro. 

Within the approximation of Eq. (2.10), the LE is expressed as a double integral containing two trajectories, 

(Po - Po) 2 + (Ps> - Po) 2 • 

(2.12) 
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As in Ref. 18, we can decompose the LE as 



M(t) = M nd (t) + M d {t) , (2.13) 

where the first term (non-diagonal) contains trajectories s and s' exploring different regions of phase space, while 
in the second (diagonal) s' remains close to s. Such a distinction is essential when considering the effect of the 
perturbation over the different contributions. 



C. Quenched disorder as a perturbation 

In order to calculate the different components to the LE (Eqs. (2.12) and (2.13)) we need to characterize the 
perturbation S. One possible choice 18 is a quenched disorder given by iVj impurities with a Gaussian potential 
characterized by the correlation length £, 



s = v(t) = ]T 



(2^ 2 ) d/2 



exp 
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The independent impurities are uniformly distributed (at positions R Q ) with density rii = A/j/V, (V is the sample 
volume). The strengths u a obey (u Q w^) = u 2 8 a p. The correlation function of the above potential is given by 



Cv(|q-q'|) = <nq)V(q')> = 



U 2 Ui 



( 47r £2)d/2 



exp 



(2.15) 



The perturbation (2.14) does not lead to the well-known physics of disordered systems, since the potential V is 
not part of Ho, but of S. Then, it acts only in the backwards propagation of the LE setup. On the other hand, 
the analogy with standard disordered systems is very useful for the analytical developments. The finite range of the 
potential allows to apply the semiclassical tool (provided ^> 1), as has been extensively used in the calculation 
of the orbital response of weak disordered quantum dots 44-47 . The finite range of the potential is a crucial ingredient 
in order to bridge the gap between the physics of disordered and dynamical systems 33,46 and to obtain the Lyapunov 
regime 18 . Moreover, taking a finite £ is not only helpful for computational or conceptual purposes, but it constitutes an 
appropriate approximation for an uncontrolled error in the reversal procedure 7io — ► — Ho+S as well as an approximate 
description for an external environment. Without entering into a discussion about what kind of perturbation more 
appropriately represents an external environment, it is reasonable to admit that the interaction with the environment 
will not be local (or short range), but will extend over certain typical length. 

Another important point to discuss concerning the appropriateness of our perturbation toward the representation of 
an external environment, is its time dependence. Taking a quenched disorder perturbation that only acts in the second 
half of the time evolution, represents a very crude approximation to the dynamics of a more realistic environment 48 . 
Moreover, since the disorder is quenched, there is no feed-back of the system on the environment. In view of our 
main result, the robustness of the Lyapunov regime respects to the details of the perturbation, this limitation should 
not prevent us of extrapolating our results to realistic cases, and envisioning their experimental consequences. In 
any case, a semiclassical approach to time dependent perturbations shows that the results of Ref. 18 remains fairly 
unchanged 49 . 
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As discussed in the previous chapter, in the leading order of h and for sufficiently weak perturbations, we can 
neglect the changes in the classical dynamics associated with the disorder. We simply modify the contributions to the 
semiclassical expansion of the LE associated with a trajectory s (or in generally to any quantity that can be expressed 
in terms of the propagators) by adding the extra phase AS of Eq. (2.11). For the perturbation (2.14) we can make 
the change of variables q = vi and write 



AS S = - — f V(q) dq . (2.16) 

The integration is now over the unperturbed trajectory C c s , and we have assumed that the velocity along the trajectory 
remains unchanged respect to its initial value vo = po/m = L s /t. 

For trajectories of length L s » £, the contributions to AS from segments separated more than £ are uncorrelated. 
The stochastic accumulation of action along the path can be therefore interpreted as determined by a random-walk 
process, resulting in a Gaussian distribution of AS S (L S ). This has also been verified numerically in Ref. 35. The 
integration over trajectories represents a for of average for exp[jrAS s ]. The ensemble average over the propagator 
(2.4) (or over independent trajectories in Eq. (2.12)) is then obtained from 



(2.17) 

and therefore entirely specified by the variance 

(AS 2 ) = -I / dq [ dq' <V(q)W)> • (2-18) 

Since the length L s of the trajectory is supposed to be much larger than £, the integral over q — q' can be taken 
from — oo to +oo, while the integral on (q + q')/2 gives a factor of L s . We thus have 

(AS 2 ) = ^ fdqC( q ), (2.19) 

resulting in 
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Where, in analogy with disordered systems 45,46 , we have defined the typical length over which the quantum phase is 
modified by the perturbation as 



dq C(q) 



u 2 rii 
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The "elastic mean free path" i and the mean free time f = £/v associated with the perturbation 50 will constitute a 
measure of the strength of the coupling. 

Taking impurity averages is technically convenient, but not crucial. Results like that of Eq. (2.20) would also arrive 
from considering a single impurity configuration and a large number of trajectories exploring different regions of phase 
space. 



D. Loschmidt echo in a classically chaotic system 

Once we have settled the hypothesis with respect to the perturbation, we can go back to Eqs. (2.12) and (2.13) 
calculate the two contributions to the Loschmidt echo. 

In the non-diagonal term the impurity average can be done independently for s and s', since the two trajectories 
explore different regions of phase space. Therefore, upon impurity average the non-diagonal term becomes 



7 



M nd (t) = \(m(t))\< 



a 



drj] C s cxp 



-jj2 (P--P0) 



(exp 



h 



AS, 



(2.22) 



We have kept the same notation for the averaged and individual LE, in order to simplify the notation, and because 
it will be demonstrated that this distinction is not crucial. According to Eq. (2.20) we have 18 
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This term depends on the perturbation, through £, and can be interpreted as a Fermi golden rule result 22 . 

In the diagonal term the trajectories s and s' of Eq. (2.12) remain close to each other. The existence of such 
types of trajectories is based on the structural stability of the manifold 27 ' 35 (opposed to the exponential sensitivity 
of individual trajectories). The actions AS, and AS S > accumulated by effect of the perturbation cannot be taken as 
uncorrelated, like in the previous case. A special treatment should be applied to the terms arising from s ~ s'. The 
small difference between s and s' is only considered through the difference of actions, and therefore 



Since s and s' are nearby trajectories, we can write 
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The difference between the intermediate points of both trajectories can be expressed using the matrix B of Eq. (2.5): 



q s (i) - q,,(t) - B-\t) (p s - p 8 ,) = B-\F)B(t) (r - r') . (2.26) 

In the chaotic case the behavior of B~ 1 (t) is dominated by the largest eigenvalue e A *. Therefore we make the 
simplification B~ 1 (i)B(t) = cxp [A(t — t)]I, where / is the unit matrix and A the mean Lyapunov exponent. Here, 
we use our hypothesis of strong chaos which excludes marginally stable regions 51 with anomalous time behavior. 
Assuming a Gaussian distribution for the random variable AS S — AS S ' 35 , in analogy with Eq. (2.17), we have 



[cxp 



-(AS S -AS S ,) 



= cxp 



~i I di l d *' ex P[ A ( < "+ < "'- 2 *)] CW(|g a (t) -«,(?)!) (r-r') 



(2.27) 



Unlike the non-diagonal case, that was obtained through the correlation potential (Eq. (2.15)), we are now led to 
consider the "force correlator" 
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Using the fact that C V y is short-ranged (in the scale of £), and working in the limit Xt » 1, the integrals of Eq. (2.27) 
yield 
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Therefore, we have 
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A Gaussian integration over (r — r') results in 
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The factor C 2 reduces to C s when we make the change of variables from r to p. In the long-time limit C s 1 oc e At , 
while for short times Cj 1 — (t/m) d . Using a form that interpolates between these two limits we have 



M d (t) 



dp ( 



27rft 2 
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with A = [am/(A 1 / 2 t)] d . Since the integral over p is concentrated around po, the exponent A is taken as the phase- 
space average value on the corresponding energy shell. The coupling E appears only in the prefactor (through A) and 
therefore its detailed description is not crucial in discussing the time dependence of M d . 

The limits of small t and weak E yield an infinite A, and thus a divergence in Eq. (2.31). However, our calculations 
are only valid in certain intervals of t and strength of the perturbation. The times considered should verify v§tjl > 1. 
Long times, resulting in the failure of our diagonal approximations (Eqs. (2.12) and (2.24)) or our assumption that the 
trajectories are unaffected by the perturbation, are excluded from our analysis. Similarly, the small values of E are 
not properly treated in the semiclassical calculation of the diagonal term M d (t), while for strong E the perturbative 
treatment of the actions is expected to break down and the trajectories are affected by the quenched disorder. This 
last condition translates into a "transport mean- free-path" 45,46 £ tT = A(k£) 2 £ much larger than the typical dimension 
R of our system. In the limit fc£ 3> 1 that we are working with, we are able to verify the condition £ tI ^> R ^> I. 

Within the above limits, our semiclassical approach made it possible to estimate the two contributions of Eq. (2.13) 
to M(t). The non-diagonal component M nd (t) will dominate in the limit of small t or E. In particular, such a 
contribution ensures that Ms =0 (t) = 1 (sec Eq. (2.9)), and that Ms(t = 0) = 1. The diagonal term will dominate 
over the non-diagonal one for perturbations strong enough to verify 



I < ^ . (2.32) 
A 

This crossover condition is extremely important, and will be discussed in detail in the sequel. 

It is worth to notice that the width a of the initial wave-packet is a prefactor of the diagonal contribution. The 
non-diagonal term, on the other hand, is independent on the initial wave-packet. Therefore, as explained in Ref. 41, 
changing our initial state (2.1) into a coherent superposition of N wave-packets would reduce M d by a factor of N 
without changing M nd . The localized character of the initial state is then a key ingredient in order to obtain the 
universal behavior. 



III. LOSCHMIDT ECHO IN THE TWO-DIMENSIONAL LORENTZ GAS 



A. Classical dynamics of Ha 



We consider in this section the case where the system Hamiltonian Ho represents a two dimensional Lorentz gas, 
i.e. a particle that moves freely (with speed v) between elastic collisions (with specular reflections) with an irregular 
array of hard disk scatterers (impurities) of radius R. Such a billiard system is a paradigm of classical dynamics, 
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FIG. 1: Schematics of a Lorentz gas showing the dispersion of two trajectories initially close to each other (with a difference 5x 
in the impact parameter x). The angle 89 between the two trajectories increases after each collision as described in the text. 

and has been proven to exhibit mixing and ergodic behavior, while its dynamics for long distances is diffusive 52-54 . 
The existence of rigorous results for the Lorentz gas has made of it a preferred playground to study the emergence 
of irreversible behavior out of the reversible laws of classical dynamics 53 . Moreover, anti-dot lattices defined in a two 
dimensional electron gas 38-40 constitute an experimentally realizable quantum system where classical features have 
been identified and measured. We will use the terms anti-dot and disk indistinctly. 

In our numerical simulations we are limited to finite systems, therefore we will work in a square billiard of area L 2 
(with N scatterers), and impose periodic boundary conditions. The concentration of disks is 



We require that each scatterer has an exclusion region R c from its border, such that the distance between the centers 
of any pair of disks is larger than a value 2R e > 2R. Such a requirement is important to avoid the trapping of the 
classical particle and the wave-function localization in the quantum case. The anti-dots density is set to be roughly 
uniform, and the concentration is chosen to be the largest one compatible with the value of i? , obtained numerically 
as c = 0.77ri? 2 /4i?2. 

The Lorentz gas has been thoroughly studied 53 , and we will not discuss here its classical dynamics in detail. We 
will simply recall some of its properties that will be used in the sequel, and present the numerical simulations that 
allow us to extract some important physical parameters. 

The chaotic character of the dynamics is a consequence of the de-focusing nature of the collisions. As illustrated in 
Fig. 1, a particle with impact parameter x will be reflected with an angle 



c = NttR 2 /L 2 . 



(3.1) 



8 = ir — 2 arctan 



x 



(3-2) 
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If we consider a second particle with impact parameter x + Sx, its outgoing angle will be 9 + 59, with 



59 = 



2 



5x. 



(3-3) 



VR 2 - x 2 



The separation between the two particles that have traveled a distance s after the collision will grow as 




(3.4) 




FIG. 2: Histogram of the distances between collisions with the disks, used in order to obtain numerically the mean free path I 
for the Lorentz gas. The solid line represents Eq. (3.5) and the dash-dotted vertical line the cut-off distance 2(R C — R). 



The next collision will further amplify the separation, due to the new impact parameters and the different incidence 
angles. 

Within the above restrictions, the exclusion distance R c completely determines the dynamical properties of the 
Lorentz gas. Among them, we are interested in the Lyapunov exponent (measuring the rate of separation of two 
nearby trajectories), the elastic mean free path I (given by the typical distance between two collisions), and the 
transport mean free path £ tT (defined as the distance over which the momentum is randomized and the dynamics can 
be taken as effectively diffusive). 

A shifted Poisson distribution 



«P[-/('-2(*.-*))] ■ lis> 2 { R e -R), 



P(s)= { (£-2(R e -R))exp[-2(R e -R)/(t-2(R e -R))] Vc " (3.5) 

if s < 2(R C - R) , 

is a reasonable guess for the distribution of lengths between successive collisions, which yields (s) = £ = v/t c , and 
is consistent with numerical simulations in the range of anti-dot concentration that we are interested in (see Fig. 2). 
Since velocity both vq and momentum po are conserved within this model we will drop their subindex. 

The elastic mean free path that we obtain from Fig. 2 compares favorably with a simple estimation of the mean 
free distance in a strip of length L and width 2R with 2cL/irR disks, 

ttR _ ttR, _ _R|_ _ nR 

2c 2 0.35i? 2 V ' ; 

The diffusive character of the Lorentz gas can be put in evidence from the time evolution of the root mean square 
displacement over a collection of initial conditions. We numerically obtain (r 2 (t)) = 2dDt (with d = 2). r tr = £t T /v 
is the mean time required to randomize the direction and D = v£ tI /2d is the diffusion coefficient. The difference 
between £ and £ tr arises from the angular dependence of the scattering cross section. Taking this factor into account 
we obtain a ratio £tr/£ which is in good agreement with the one obtained from the independently determined £ and 

£tr- 

There are known various estimations of the Lyapunov exponent of the Lorentz gas in different regimes. Considering 
the three-disk problem, Gaspard et al. 55 obtained 



A = 



2R V - 2R 



In 
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FIG. 3: Lyapunov exponent A of the Lorentz gas as a function of the mean free path I. The black dots represent our numerical 
values and the solid line the analytical estimate of Eq. (3.12). The dashed line indicates Laughlin's approximation (Eq. (3.8) 
and the open dots are the quantum values obtained from the decay of the LE (Fig. 6 in the sequel). Inset: the same plot in 
log-log scale highlighting the agreement between the different approximations in the region of very small concentrations (large 
I). 



Considering a periodic Lorentz gas (repeated Sinai billiard) Laughlin proposed the form 10 



(3i 



where (5 is a geometrical factor of order 1. In the diluted limit (c <C 1), van Beijeren and Dorfman 56 showed that 



N 

A = 2 Rv ( I - hi 2 - 0..177 - In 



NR 2 



L 2 



(3.9) 



Numerically, the procedure of Benettin et al. 57 is usually followed in order to obtain Lyapunov exponents. Two 
nearby trajectories are followed, and their separation is scaled down to the initial value Sxo after a characteristic time 
t (that we take it to be larger than the collision time). We can then limit the numerical errors, and avoid entering the 
diffusive regime where two initially close trajectories follow completely independent paths. The Lyapunov exponent 
results from the average over the expanding rates in the different intervals, 



A = lim V - Vl 

n^oo n ~ — ' Si 



In 



Sx 



(3.10) 



where Sj is the length of the j-th interval, and Sxj the separation just before the normalization. Technically, we 
should work with distances in phase-space, rather than in configuration space, but the local instability makes this 
precision unnecessary. 

Benettin's algorithm can also be used for a semi- analytical calculation of the Lyapunov exponent. Taking the length 
distribution of Eq. (3.5) to obtain the average separation after a collision from Eq. (3.4), and identifying the average 
over pieces of the trajectory with a geometrical average over impact parameters, we can write 



A = 



dx In 



1 + 
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Performing the integration yields 
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As shown in Fig. 3, the above expression reproduces remarkably well the numerical calculations of the Lyapunov 
exponent. It agrees also with the result of van Bcijcrcn and Dorfman in the dilute limit, and gives good agreement 
with Laughlin's estimation. 



B. Perturbation Hamiltonian 



The quantum and classical fidelity measure the sensitivity of a given system to a perturbation of its Hamiltonian. In 
Ref. 18 a quenched disorder environment was taken as the perturbation and the relaxation rate was found to depend 
only on the system Hamiltonian; the details of the perturbation are not important beyond some critical strength. 
Subsequent works have tried this perturbation 25 and others 20-33 confirming the universality of the result of Ref. 18. 
It is also useful to verify such an universality by considering completely different perturbations on a given system. 
The Lorentz gas is an ideal case, since it can be perturbed by a quenched disorder or by the distortion of the mass 
tensor, introduced in Ref. 21 and briefly discussed in the sequel. 

The isotropic mass tensor of 7i , of diagonal components m , can be distorted by introducing an anisotropy such 
that m xx = mo(l + a) and m yy = mo/(l + a). This perturbation is inspired by the effect of a slight rotation of the 
sample in the problem of dipolar spin dynamics 58 , which modifies the mass of the spin wave excitations. The kinetic 
part of the Hamiltonian is now affected by the perturbation, that writes as 

^ a) = a A 2_ JL. (3.i3) 

v ; 2m 1 + a 2m 
In our analytical work we will stay within the leading order perturbation in a. That is, 

= ^ (Pl-Pl)- (3.14) 

Making the particle "heavier" in the x direction (i.e. we consider a positive a) modifies the equations of motion 
without changing the potential part of the Hamiltonian. It is important to notice that, unlike the case of quenched 
disorder, the perturbation (3.13) is non-random, and will not be able to provide any averaging procedure by itself, 
but through the underlying chaotic dynamics. 

Numerical simulations of the evolution of two trajectories with the same initial conditions, the first one governed 
by Tlo and the second one by Ho + E, show that the distance in phase space grows exponentially with the same 
Lyapunov exponent that amplifies initial distances. The classical dynamics is then equally sensitive to changes in the 
Hamiltonian as to changes in the initial conditions 59 . 

For a hard wall model, like the one we are considering, the perturbation (3.13) is equivalent to having non-specular 
reflections. One can resort to the minimum-action principle (see appendix A) to obtain a generalized reflection law: 



v x (m x nl - m y n 2 x ) - 2v y m y n x n y 

v = 

m x n y + m y n x 

v y (m y nl - m x nl) - 2v x m x n x n y 

v = - 

y m x n y + m y n x 



(3.15a) 
(3.15b) 



where v' x and v' are the two components of the velocity after a collision against a surface defined by its normal 
unitary vector (n x ,n y ). Eqs. (3.15) allow to show that the distortion of the mass tensor is equivalent to an area 
conserving deformation of the boundaries as x —> x(l + £), y — > y/(l + £), as used in other works on the LE 24 , where 
£ = VI + a — 1 is the stretching parameter. An illustrative example of the equivalence of the perturbations is shown 
for a stretched stadium billiard in Fig. 4. In panel a, a typical trajectory is shown; while in panels b and c we see the 
trajectories resulting of a perturbation of the mass tensor or the stretch of the boundaries respectively. 



C. Semiclassical Loschmidt echo 



We calculate in this chapter the Loschmidt echo for the system whose classical counterpart was previously discussed; 
Tto describes a Lorentz gas and £ is given by Eq. (3.13). We adapt to the present perturbation the semiclassical 
approach of Sect. II and Ref. 18. As before, we take as initial state a Gaussian wave-packet of width a (Eq. (2.1)). 

The semiclassical approach to the LE under a weak perturbation S is given by Eq. (2.12), with the extra phase 
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FIG. 4: A single typical trajectory on the (hard wall) Bunimovich stadium, the arrows indicate the initial and final positions of 
the particle, (a) Normal stadium (b) Stretched stadium, (c) Normal stadium and particle with a perturbed mass, Eqs. (3.13). 
The strengths of the stretching and the mass distortion are such that the initial condition renders the same trajectory in both 
cases, illustrating the corresponding map (see text). 



AS S = f dt £ a (q(t),q(t)). (3.16) 
Jo 

The sign difference with Eq. (2.11) is because the perturbation is now in the kinetic part of the Hamiltonian. On the 
other hand, as explained before, this sign turns out to be irrelevant. 

With the perturbation of Eq. (3.14) we have to integrate a piecewise constant function (in between collisions with 
the scatterers), obtaining 

AS S = ^ f>(2<-, 2 ). (3.17) 

i=l 

We have used v\ + v y = v 2 , and have defined n as the free flight time finishing with the i-th collision, v Vi is the y 
component of the velocity in such an interval, and N s as the number of collisions that suffers the trajectory s during 
the time t. 

As we saw in chapter III A, the free flight times Ti (or the inter-collision length vTi) have a shifted Poisson distribution 
(Eq. (3.5) and Fig. 2). This observation will turn out to be important in the analytical calculations that follow since 
the sum of Eq. (3.17) for a long trajectory can be taken as composed of uncorrelated random variables following 
the above mentioned distribution. It is important to remark that, unlike the case of Sec. II, the randomness is not 
associated with the perturbation (which is fixed), but with the diffusive dynamics generated by Tin- 



D. Non-diagonal contribution 



As in the case of Sec. II, the non-diagonal contribution is given by the second moment 

2 2 / N * \ 



(3.18) 



Separating in diagonal (i = j) and non-diagonal (i ^ j) contributions (in pieces of trajectory) we have 



(AS 2 ) 



(t?) (4 «> - 4, 2 (vl) + , 4 ) + (N s 1) (nf (4 «> 2 - 4, 2 «) + , 4 ) 



(3.19) 



We have assumed that different pieces of the trajectory (i ^ j) are uncorrelated, and that within a given piece i, Ti 
and v Vi are also uncorrelated. According to the distribution of time-of-flights (3.5) we have 
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(r) 
(r 2 ) = 2r ( 



2 



(3.20a) 
(3.20b) 



Assuming that the velocity distribution is isotropic (P(9) = l/2n, where 9 is the angle of the velocity with respect 
to a fixed axis) is in good agreement with our numerical simulations, and results in 



<, 2 )=, 2 <sin 2 tf) = - 
«>=, 4 <sin 4 0> = 



v" 
2 

3« 4 



(3.21a) 
(3.21b) 



We thus obtain that 4 (v 2 .) 2 — 4v 2 (v 2 .) + v A = 0, implying a cancellation of the cross terms of (A5 2 ), consistently 
with the lack of correlations between different pieces that we have assumed. We therefore have 



^ AS 2^ = cPmlNsT'Zv 4 ' 



(3.22) 



For a given t, N s is also a random variable, but for f > r c we can approximate it by its mean value t/r c and write 

(A*?) = ■ (3.23) 

We therefore have for the average echo amplitude 



(m(t)} ~ exp 



2 2 4 _t 



m 2 



2 \d/2 

irf? ) 



exp 
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vt 


= exp 




~2~l 



(3.24) 



where we have again used C s as a Jacobian of the transformation from r to p s and we have defined an effective mean 
free path of the perturbation by 



mlv 2 £ 2 



(3.25) 



The effective mean free path £ = v f should be distinguished from I = vt c since the former is associated to the 
dynamics of S and Tlo, while the latter is only fixed by TCq. Obviously, our results are only applicable in the case of 
a weak perturbation verifying I » i. From Eq. (3.24) we obtain the non-diagonal component of the LE as 



M^{t) = \( m {t))Y 



exp 



vt 

T 



(3.26) 



In the next chapters we study the conditions under which the correlations not contained in the FGR approximation 
dominate the LE, while in Sec. IV we will test the above results against numerical simulations. 

E. Diagonal contribution 

As in Sec. II, we have to discuss separately the contribution to the LE (Eq. (2.12)) originated by pairs of trajectories 
s and s' that remain close to each other. In that case the terms AS S and AS S > are not uncorrelated. The corresponding 
diagonal contribution to the LE is given by Eq. (2.24), and then we have to calculate the extra actions for s ~ s' . As 
in Fig. 1, we represent by 9 (9 + S) the angle of the trajectory s (s 1 ) with a fixed direction (i.e. that of the x-axis). 
We can then write the perturbation (Eq. (3.13)) for each trajectory as 



S, = 



2m 
a 

2m 



p 2 (2 sin 2 9 - 1) , 

p 2 (2 sin 2 9 - 25 sin 29 - 1) + 0{S 2 ) 



(3.27a) 
(3.27b) 
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Assuming that the time-of-flight n is the same for s and s' we have 

A*. - A,S'. S - f dt 5(t) sin [20(t)] 



wo Jo 



(3.28) 



The angles (5 alternate in sign, but the exponential divergence between nearby trajectories allows to approximate 
the angle difference after n collisions as \S n \ = |(5i|e AnT °. A detailed analysis of the classical dynamics shows that 
the distance between the two trajectories grows with the number of collisions as c?i = \Si\ vt\, d 2 — d\ + |<5 2 | vt 2 , and 
therefore 
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By eliminating \5\\ we can express an intermediate angle 5(t) as a function of the final separation |r — r'| = d^ s 



8(t) 



|r-r'| e AT <= - 1 



XL 
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1 



(3.30) 



where again we have used that t = N s t c is valid on average. Assuming that the action difference is a Gaussian random 
variable, in the evaluation of Eq. (2.24) we only need its second moment 



(AS S — AS S >) ) ~ a rriQV 



r . -2,4 |r-r'| 2 ( e Are - V 1 
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^) (Jo ^Jo ^ ' e>J+>1 ' sin W)] sin [20(7')]^ . (3.31) 



As before, we assume that the different pieces are uncorrelated and the angles 0i uniformly distributed. Therefore 
(sin [20*] sin [26 j] > = % /2 and 
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where we have taken the limit Xt >• 1, and defined 
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(3.32) 
(3.33) 
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2 Xt ) e A ^ + l 



(3.34) 



Our result (3.33) is analogous to Eq. (2.29) obtained in the case of a perturbation by a quenched disorder. Obviously, 
the factor A is different in both cases, but we use the same notation to stress the similar role as just a prefactor of 
M d . Performing again a Gaussian integral of M d over r — r' we obtain 



2irh 2 



exp 



2o 2 2 
-p- (P« -Po) 



(3.35) 



Under the same assumptions than in Sec. II D, we are lead to a result equivalent to that of Eq. (2.31), 



M d {t) ~ Ae- Xt , (3.36) 

with A = [amo/(A 1 ^ 2 t)] d . Therefore, for long times the diagonal part of the Loschmidt echo decays with a rate given 
by the classical Lyapunov exponent of the system, 



lhn(-i ln[M^)])=A. 



(3.37) 



Of course this limit actually means t 3> 1/A but still lower than the time at which cither localization or finite size effect 
appears. In the next chapter we will study the competition between the diagonal and non-diagonal contributions. 
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F. Diagonal vs. non-diagonal contributions 

As we have previously shown, the Loschmidt echo is made out of non-diagonal and diagonal components, and within 
the time scales above specified, it can be written as 



+ .4 exp [-At] . (3.38) 

Such a result holds for the perturbation S that we have discussed in this section (Eq. (3.13)), as well as for the 
quenched disorder of Sec. II (Eq. (2.14). The only difference lays in the form of the "elastic mean free path" I and 
the prefactor A, both of which are perturbation dependent. The decay of the LE will be controlled by the slowest of 
the two rates. A weak perturbation implies £ > v/X and a dominance of the non-diagonal term, while for sufficiently 
strong perturbations verifying I < v/X (but weak enough in order not to modify appreciably the classical trajectories), 
the diagonal term (governed by the Lyapunov exponent) sets the decay of the LE. In Ref. 22 the regime of dominance 
of the non-diagonal and diagonal component has been respectively interpreted and referred to as a Fermi Golden Rule 
Lyapunov regimes and we will use both terminologies in the discussions that follow. 

The Lyapunov regime is remarkable in the sense that its decay rate is an intrinsic property of the system and does 
not depend on the perturbation that gives rise to the decay. This behavior, predicted in Ref. 18 has been observed in 
numerical simulations done on a number of systems 20 ~ 33 . 

From the previous discussion it is clear that the Lyapunov regime can only be observed beyond a critical value 
of the perturbation. The condition stated above for the strength of the perturbation, along with Eq. (3.25), yields 
for the model discussed in this section a critical value of the perturbation parameter a beyond which the Lyapunov 
regime is obtained, 



M(t) = exp 



vt 

T 



a c = 




to 



(3.39) 



We will discuss in Sec. IV the physical consequences of the above critical value and its dependence on various physical 
parameters. 

We finish this chapter with the discussion of the perturbation dependent (non-diagonal) regime. In this Fermi 
Golden Rule regime the LE is equal to the return probability P(t) = |(f/'o| exp [— i (7Y + S) t/H] \tpo}\ 2 , whose decay 
rate does not show saturation at the Lyapunov exponent but rather follows Eq. (3.25) for the whole range of 
parameters 60 . The full connection between the exponent of Eq. (3.25) and the one we would get from a complete 
FGR approach, was clarified using a random matrix treatment 25 . In fact, a rough estimation of the Fermi golden 
rule is obtained considering a particle moving along the principal axis to* of a square box of sides L x = L y = I. The 
available density of states ^ , corresponds to a 1-d tube of length i. Hence 

1 = 2tt 2 
T ft 1 1 A 



2tt 

Y 



p- 



2m 



£m 
it hp 



in agreement with the semiclassical calculation of Eq. (3.25). Of course, this estimation does not make justice to 
the chaotic nature of Ho- This breaks the selection rules of our simple perturbation and enables a random matrix 
approximation for E, that mixes eigenstates of Ho which follow a Wigner-Dyson statistics. Hence the perturbation 
breakdown of the FGR regime can not be calculated from the parameters introduced above. 

The LE has also been studied in systems where the perturbation structure prevents the application of the Fermi 
Golden Rule. The result is that in general the decay rate of the LE before the Lyapunov regime is given by the width 
of the local density of states of the perturbation, which for particular systems coincides with the exponent given by 
the FGR 24 . 
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FIG. 5: Time decay of the Loschmidt echo M(t) for different values of the perturbation strength (a) and concentration of 
impurities (c). (a) c = 0.157 and a =0.004, 0.007, 0.01, 0.015, 0.02, 0.03, 0.05, 0.07, 0.1 (from top to bottom); (b) c = 0.195 
and a =0.004, 0.007, 0.01, 0.015, 0.02, 0.03, 0.04, 0.05, 0.06, 0.07, 0.08, 0.1, 0.15; (c) c = 0.289 and a =0.004, 0.007, 0.01, 
0.015, 0.02, 0.03, 0.04, 0.05, 0.06, 0.07. The time is measured in units of h/V, where V is the hopping term of the tight-binding 
model (see appendix B). The doted lines represent the best fits to the decay, as described in the text. 

IV. UNIVERSALITY OF THE LYAPUNOV REGIME 
A. Correspondence between semiclassical and numerical calculations 

The semiclassical results obtained in the previous sections are valid in the small wavelength limit, and relay on 
various uncontrolled approximations. It is then important to perform numerical calculations for our model system in 
order to compare against the semiclassical predictions, and to explore parameter regimes inaccessible to the theory. 
In this section we use the same numerical method of Ref. 21 for the Lorentz gas (described in detail in appendix B), 
and extend the results in order to sustain the discussion on the universality of the Lyapunov regime. We will first 
focus on the behavior of the ensemble averaged Loschmidt echo, followed by a thorough discussion of the averaging 
procedure and the individual behavior. 

We typically worked with disks of radius R = 20a , and with a Fermi wavelength Xp = 2tt /kp = 16/3a. Here, 
a is the irrelevant lattice unit of our tight-binding model, which is decreased until the results only depend on the 
relation between physical parameters. The smallest system-size allowing to observe the exponential decay of M(t) 
over a large interval was found to be L = 200a which means the consideration of a Hilbert space with 4xl0 4 states. 
We calculated M(t) for different strengths of the perturbation a and concentration of disks c. In Fig. 5 we show our 
results for c = 0.157, 0.195 and 0.289 (panels a to c respectively), and different values of a. 

The time evolution of the LE presents various regimes. Firstly, for very short times, M(t) exhibits a Gaussian decay, 
M(t) = cxp [— ba 2 t 2 ] , where b is a parameter that depends on the initial state, the dynamics of Ho and the form of 
the perturbation S. This initial decay corresponds to the overlap of the perturbed and unperturbed wave-packets 
whose centers separate linearly with time by the sole effect of the perturbation. This regime ends approximately at 
the typical time of the first collision. 

Secondly, for intermediate times we find the region of interest for the semiclassical theory. In this time scale the 
LE decays exponentially with a characteristic time r^. We reserve the symbol for the decay rate, in view of its 
interpretation in terms of quantum decoherence (as we discuss in Sec. V). For small perturbations, depends on a. 
We observe that for all concentrations there is a critical value a c beyond which is independent of the perturbation. 
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FIG. 6: Extracted values of the decay rate of the LE as a function of the perturbation strength a for the three concentrations 
of Fig. 5. The rates are normalized to the group velocity of the initial wave-packet v (1/(vt$)) is given in units of a -1 ; c = 0.157 
(circles), 0.195 (squares) and c = 0.289 (triangles). The solid lines are the corresponding classical Lyapunov exponents and the 
dashed lines are fits to the quadratic behavior predicted by Eq. (3.25). The predicted coefficients for the three concentrations 
are 72a - , 55a -1 and 33a -1 , while the obtained ones are 92a~ , 50a~ and 37a~ respectively. In the inset, a log-log scale of 
the same data to show the quadratic increase of l/r^ for small perturbations. 



Clearly, the initial perturbation-dependent Gaussian decay prevents the curves to be superimposed. 

Finally, for very large times the LE saturates at a value that depends on the system size L. This regime 
is discussed in detail in the next chapters. However, let us observe that in the crossover between the exponential 
decay and the long time saturation there is a power-law decay with a perturbation independent exponent. This is 
a manifestation of the underlying diffusive dynamics that leads to the isotropic state. Therefore, it could be related 
to the Ruelle-Perricot resonances of the classical Perron-Frobenius evolution operator used to calculate a classical 
version of the LE 32 . 

In order to compare our numerical results of M(t) with the semiclassical predictions, we extract by fitting In M(t) 
to In [A exp(— t/r^) + M^] . The dashed lines in Fig. 5 correspond to the best fits obtained with this procedure. The 
values of for the different concentrations are shown as a function of the perturbation strength in Fig. 6. In agreement 
with our analytical results of the previous section, we see that 1/t^ grows quadratically with the perturbation strength 
up to a critical value a c , beyond which a plateau appears at the corresponding Lyapunov exponent. The dashed 
lines are the best fit to a quadratic behavior. The values obtained in this way agree with those predicted by the 
semiclassical theory (Eq. (3.25)) for the non-diagonal (FGR) term. The saturation values above a c are well described 
by the corresponding Lyapunov exponents (solid lines), in agreement with the semiclassical prediction (Eq. (3.37)). 
The very good quantitative agreement between the semiclassical and numerical calculations for the Lorentz gas (as 
well as in the case of other models 22,24 ' 25 ) strongly supports the generality of the saturation of at a critical value 
of the perturbation strength. 

The FGR exponent, which depends on H.q but not much on its chaoticity 30 , is given by the typical squared matrix 
element of S, and the density of connected final states 1/A. Hence, different Ho change the wave-functions. That is 
why we observe that, for fixed perturbation strength a, the factor v/£ depends on the concentration of impurities of 
Tio (see inset of Fig. 6, where a log- log scale has been chosen in order to magnify the small perturbation region). 

Notably, the dependence of vfl with 7io leads to a counter-intuitive effect (clearly observed in the inset of Fig. 6), 
namely that the critical value needed for the saturation of 1/r^ is smaller for less chaotic systems (smaller A). The 
reason for this is that in more dilute systems E is constant over larger straight pieces of trajectories (in between 
collisions) , leading to a larger perturbation of the quantum phase and resulting in a stronger effective perturbation. 



B. Universality of the Lyapunov regime in the semiclassical limit 



Our semiclassical analysis yielded a critical value of the perturbation to enter in the Lyapunov regime (Eq. (3.39)), 
that vanishes in the semiclassical limit, a c — > for h (or Xp) — > 0, implying the collapse of the Fermi Golden Rule 
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FIG. 7: Decay rates 1/t> for different wavelengths Af of the initial wave-packet for a concentration c = 0.195 with the same 
units as in Fig 6. Solid line: classical Lyapunov exponent. Dashed line: the FGR quadratic behavior. Note that for decreasing 
Af the critical perturbation diminishes, implying a collapse of the Fermi Golden Rule regime. 



regime. This behavior is reproduced by our numerical calculations (Fig. 7). There, we decreased Af while keeping 
fixed the size a of the initial wave packet. A point that should not be over-sighted is that the perturbation S (Eq. 
(3.13)), for a given value of the parameter a, scales with the energy in a way that the underlying classical trajectories 
are always affected in the same way by the perturbation. The extracted crossover values of a c are in quantitative 
agreement with Eq. (3.39), decreasing with Af in the interval that we were able to test. 

Other choices of the perturbation S, such as the quenched disorder of Refs. 18 and 25, can be shown to give 
critical values that decrease with decreasing h as in Eq. (3.39), provided that the perturbation is scaled to the proper 
semiclassical limit. That is, for a fixed perturbation potential, we should take the limit of Xp — > 0. As a result, if 
we keep H constant and decrease Af by increasing the particle energy, we should scale up the perturbation potential 
consistently (assuming that TLq generates the same dynamics at all energies). 

We conclude that, in the semiclassical limit, any perturbation will be strong enough to put us in the Lyapunov 
regime, in consistency with the hypersensitivity expected for a classical system. This is not unexpected as in this 
limit the Ehrenfest time diverges and the correspondence principle should prevail. 

We can draw a the critical perturbation strength separating FGR from Lyapunov regimes versus a scaling parameter 
determined by the particle energy (or inverse h), as shown in Fig. 8. The shaded region corresponds to the Fermi 
Golden Rule regime and the clear one to the Lyapunov regime. The line that divides both phases is given by Eq. (3.39), 
and the dots correspond to numerical values of a c extracted from Figs. 6 and 7. Of course, there is another transition 
from FGR to perturbation appearing when S ~ A, which we avoid drawing since a a *C a c . This perturbative value 
also goes to zero in the semiclassical limit of Af — > 0. Also the Lyapunov regime is bounded from above by an h 
independent critical value a p marking the classical breakdown that we discuss bellow. 

The interesting conceptual feature highlighted by Fig. 8, is the importance of the order in which we take the 
limits of £ and Af going to zero. Two distinct results are obtained for the different order in which we can take this 
double limit. As depicted in the figure (with arrows representing the limits), lim,x F _>o lim£_ >o 1/ T = 0- On the other 
hand, taking the inverse (more physical) ordering lims^o liniA F ^o = A the semiclassical result is obtained. The 
resulting "phase diagram" representation for the different regimes of the LE serves us to remember that most often 
one is working in the thermodynamic side corresponding to the Lyapunov region. 

Our semiclassical theory clearly fails when the perturbation is strong enough (or the times long enough) to appre- 
ciably modify the classical trajectories. This would give an upper limit (in perturbation strength) for the results of 
Sec. III. A more stringent limitation comes from the finite value of h, due to the limitations of the diagonal approx- 
imations and linear expansions of the action that we have relied on. In other systems, like the quenched disorder in 
a smooth stadium 25 , the upper critical value of the perturbation (for exiting the Lyapunov regime) can be related 
to the transport mean free path of the perturbation £ tT , which is defined as the length scale over which the classical 
trajectories are affected by the disorder 45 . 

We can obtain in our system an estimate of £t T by considering the effect of the perturbation on a single scattering 
event. The difference 59 between the perturbed and unperturbed exit angles after the collision can be obtained using 
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FIG. 8: Regime diagram for the Loschmidt echo as a function of the perturbation and the energy (or inverse K). The grayed 
area is the FGR regime, while the clear one is the Lyapunov regime. The line that divides both regimes is Eq. (3.39). The dots 
are the numerical values obtained from Figs. 6 and 7. The arrows schematize the possible ordering of the classical double limit 
of the perturbation and the wavelength going to zero. Notice how the lower one gives always zero while the upper (correct) 
one gives A since it remains always in the Lyapunov regime. 



Eqs. (3.15), which results in 

/ v • n\ 2 

56 ~ An x n y y J a, (4.1) 

where v is the initial velocity of the particle and n is the normal to the surface. 

Assuming that the movement of the particle is not affected by chaos (non-dispersive collisions), one can do a random 
walk approach and estimate the mean square distance after a time r tr from the fluctuations of the angle in Eq. (4.1). 
We estimate the transport mean free time as that at which the fluctuations are of the order of R, and obtain 

AR 2 

f » = s«- < 42 » 

assuming a uniform probability for the angle of the velocity as before. Eq. (4.2) is used to get the upper bound 
perturbation a p for the end of the Lyapunov plateau, 



We obtain a p ~ 0.23, 0.29 and 0.43 respectively for increasing magnitude of the three concentrations shown in Fig. 6. 
It is rather difficult to reach numerically these perturbations in our system, since the initial Gaussian decay drives 
M{t) very quickly towards its saturation value, preventing the observation of an exponential regime. Despite this 
difficulty, we observe in Fig. 7 that Lyapunov regime is a plateau that ends up for sufficiently strong perturbations. 
For the range we could explore the limiting values are in qualitative agreement with the estimation from Eq. (4.3). 



C. Ehrenfest time and thermodynamic limit 



We studied in the previous chapter the behavior of the Lyapunov regime in the semiclassical limit ft — > 0; let us 
now turn our attention to the consequences of having a finite value of ft. In this case, one expects the propagation of a 
quantum wave-packet to be described by the classical equations of motion up to the Ehrenfest time tg, after which the 
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FIG. 9: M(t) for two systems with the same number of states N but different sizes L, showing that the saturation time does 
not depend on N but on the ratio L/a. The Ehrenfest time is marked with an arrow. 



quantum-classical correspondence breaks down . Typically, tE is the time when interference effects become relevant, 
and in a classically chaotic system it typically scales as ln[7i]. 

In other systems where the Lyapunov regime of the LE has been observed, such as chaotic maps or kicked systems, 
tE coincides with the saturation time t 8 = 1/A ln[iV]. This is because in these systems the number of states N plays 
the role of an effective Planck's constant h e s — l/N. Therefore, when in these systems the LE is governed by a 
classical quantity, the whole range of interest occurs before the Ehrenfest time. This observation lead Benenti and 
Casati 26 to propose that the independence of the decay rate on the perturbation strength is a trivial consequence of 
the quantum-classical correspondence before tE- 

In the Lorentz gas, however, we can differentiate between the time scales t s and tE by appropriately controlling the 
parameters. The saturation time is given by 

ta^ln-. (4.4) 
A a 

The Ehrenfest time, defined as the time it takes for a minimal wave-packet of wavelength A^ to spread over a distance 
of the order of i? 54 , is given by 

t B ~ T In— . (4.5) 

A Ap 

Our numerical calculations support these approximations. We show in Fig. 9 M{t) for two systems with the same 
number of states N (hence same H e g), but two different sizes L. N is controlled by the discretization step a — L/N. 
The different saturation values (and times) observed imply that for the Lorentz gas t e and t s are independent of each 
other. Clearly we can study the LE for times arbitrarily larger than tp by increasing the system size L, which controls 
t s , and keeping all other parameters (including tp, marked with a dashed line in Fig. 9) fixed. We can see that the 
exponential decay of the LE continues for times larger than tp , up to the saturation time. 

The above results are further evidence of the universality of the Lyapunov regime, for it persists for arbitrarily large 
times in the thermodynamic limit of the size of the system going to infinity. This is exemplified in Fig. 10, where 
we show M(t) for increasing sizes [L — 200a, 400a and 800a) for a fixed concentration (c = 0.195) and perturbation 
(a = 0.024). Notice how the exponential decay regime extends, as L grows, for times larger than tE, where the 
correspondence principle does not prevail. The survival of a classical signature of the quantum dynamics after the 
Ehrenfest time is due to a more complex effect, namely the decoherence that washes out terms of quantum nature. 
We will discuss this process in detail in the next chapter. 

In the inset of Fig. 10 we see the saturation value M x as a function of the inverse system size 1/L 2 . This dependence 
was expected from earlier works on the LE 16 . Supposing that for long times the chaotic nature of the system will 
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FIG. 10: M(t) for different system sizes L — 200,400 and 800 showing the longer exponential regime for larger L. Inset: 
saturation value M>o as a function of inverse system size. The straight line is the fit Moo = (0.6)(a/L) 2 . 

equally mix the TV = {L/cr) 2 levels appreciably represented in the initial state with random phases <f>j, we write 

2 



Moo = "i- 

N 2 



1 

N 



(4.6) 



We also show in a straight line the best fit to the data, M m = (0.6 ± 0.1) {u/L) which confirms the prediction. 



D. Individual vs. ensemble-average behavior 



In order to make analytical progress, in our semiclassical calculations and in those of Rcf. 18, an ensemble average 
was introduced (over realizations of the quenched disordered perturbation or over initial conditions). This approx- 
imation raises the question of whether the exponential decay of M(t) is already present in individual realizations 
or, on the contrary, the averaging procedure is a crucial ingredient in obtaining a relaxation rate independent of the 
perturbation 61 . 

As it was discussed in Sects. II and III, for trajectories longer than the correlation length £ of the perturbation, 
the contributions to AS from segments separated by more than £ are uncorrelated. This leads us to consider that the 
decay observed for a single initial condition will be equivalent to that of the average. In this section we test this idea 
numerically. 

For large enough systems presenting a large saturation time, we expect M(t) to fluctuate around an exponential 
decay. This expectation is clearly supported by our numerical results shown in Fig. 11, where we present M(t) for 
three different initial conditions in a system with L = 800a and fixed a — 0.024. An exponential decay with the 
semiclassical exponent is shown for comparison (thin solid line). 

In order to obtain the exponent of the decay with a good precision, we can calculate M{t) for a single initial 
condition in a large enough system. Alternatively, our results show that it is correct to obtain the exponent through 
an ensemble average to reduce the size of the fluctuations. However, as the former method is computationally much 
more expensive, we resort to the latter. 

This situation is analogous to the classical case where one obtains the Lyapunov exponent from a single trajectory 
taking the limit of the initial distance going to zero and the time going to infinity, or else resorts to more practical 
methods like Benettin et al algorithm that average distances over short evolutions. 

Notice that in the Lorentz gas the average over initial conditions and the average over realizations of the impurities 
positions are equivalent. In all cases we have implemented the last choice for being computationally convenient, and 
we use the term initial conditions to refer also to realizations of Hq. 
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FIG. 11: M(t) for three different single initial conditions of the wave-packet. All the curves oscillate around the straight line, 
which is the decay corresponding to the Lyapunov exponent. 

In particular for our calculations, the average is constrained to those systems where the classical trajectory of the 
wave-packet collides with at least one of the scatterers. This restriction helps avoiding those configurations where a 
"corridor" exists, in which case M(t) presents a power- law decay possibly related to the behavior found in intcgrable 
systems 30 . 

E. Effect of the average procedure 

The averaging of quantities that fluctuate around an exponential decay is a delicate matter, since different procedures 
can lead to quite different results. In particular, for the LE it has been noted that averaging M(t) over initial conditions 
can result in an exponential decay different than the one of a single initial condition 23 ' 61 . This effect can be attained, 
for instance, if we suppose that M(t) for single conditions decays exponentially with a fluctuating exponent A + 5X, 
where SX is randomly distributed with uniform probability between — <t\ and a\. Given the exponential dependence 
of M(t) in A, the phase space fluctuations of the Lyapunov exponent will induce a difference between the average 
In M(t) and that of M(t). The former procedure is more appropriate in order to have averages of the order of the 
typical values. On the other hand, if the fluctuations of the exponent are small, both procedures give similar results. 
This is the situation we found in our model system. 

For the Lorentz gas we calculated (M(t)) and (lnM(t)) and extracted the decay rates of the exponential regime 
using the fit described in Sec. IV A. Typical results are shown in Fig. f 2 as a function of the perturbation. We observe 
that both averaging procedures give values of that are indistinguishable from each other within the statistical error. 

A typical set of curves of M(t) averaged following the two procedures is shown in the inset of Fig. f 2. We can 
see that, even though for the short and long time regimes the two procedures give slightly different results, the 
intermediate exponential decay regime has approximately the same decay rate in both cases. 

V. ANALYSIS OF DECOHERENCE THROUGH THE LOSCHMIDT ECHO 

A. Classical evolution of the Wigner function 

As discussed in the introduction, the Loschmidt echo can be obtained from the evolution of the Wigner function 
with the perturbed and unperturbed Hamiltonians (Eq. (1.2)). Such a framework is particularly useful in the study of 
decoherence, as the Wigner function is a privileged tool to understand the connection between quantum and classical 
dynamics 36,43 . 

The evolution of the wave-functions in terms of the propagators (Eq. (2.2)) can be used to express the time- 
dependence of the Wigner function as 
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FIG. 12: Values ol extracted from averaging the echo (M) (triangles) and the logarithm of the echo (log M) (circles) for 
c = 0.289 and L — 200, as a function of the perturbation strength. Both averaging procedures yield approximate the same 
decay rate. Inset: typical curves of (M) and (logM) for the same set of parameters. The straight dotted line represents an 
exponential decay with a rate given by the semiclassical prediction. 
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The semiclassical expansion of the propagators (Eq. (2.4)) leads to the propagation of the Wigner function by 
"chords" 62-64 , where pairs of trajectories (s, s') traveling from (r — Sf/2, r + Sf/2) to (r — Sr/2, r + Sr/2) have to be 
considered (Fig. 13). In the leading order in h we can approximate the above propagators by sums over trajectories 
going from r to r 
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where p s (p s ) and p s / (p s >) are the initial (final) momenta of the trajectories s and s', respectively. The semiclassical 
evolution of the Wigner function is given by 



W(r,p;t) = (2wh) d j dr J dp IF(r, p; 0) ]T S (p 



The dominant contribution arises from the diagonal term s = s' 



W c (r 



,p,t) = J dr J2 C s S(p-p s ) VF(r,p s ;0) 



(5.3) 



s(r,r,t) 



25 




FIG. 13: Schematics of the classical trajectories involved in the semiclassical approximation to the propagation of a Wigner 
function. 

Using the fact that C s is the Jacobian of the transformation from r to p s , we have 



(r, p; t) = J dp s S(p- p s ) W(r, p s ; 0) , 



(5.4) 



where the trajectories considered now are those that arrive to r with momentum p. We note (r,p) the pre-image of 
(r,p) by the equations of motion acting on a time t. That is, (r, p) = X t (r,p). The momentum integral is trivial, 
and we obtain the obvious result 



W c (r,p;t) = W(r,p;0) 



(5.5) 



with (r, p) = X t 1 (r, p). Since X t conserves the volume in phase-space, at the classical level the Wigner function 
evolves by simply following the classical flow. 

B. Fine structure of the Wigner function and non-classical contributions to the Loschmidt echo 

As indicated in Eq. (1.2), the Loschmidt echo is given by the phase-space trace of two Wigner functions associated 
with slightly different Hamiltonians (Ho and Ho + S). In order to facilitate the discussion, we introduce the density 
(or partial trace) /s writing the LE as 



with 
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The semiclassical evolution of /s is given by four trajectories, as illustrated in Fig. 14. 

As we have consistently done in this work, we take Gaussian wave-packet (of width a) as initial state. Its associated 
Wigner function reads 
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FIG. 14: Four classical trajectories used to compute semiclassically the Loschmidt echo through the evolution of two Wigner 
functions associated with different Hamiltonians. 



Assuming that £ constitutes a small perturbation, after a few trivial integrations we obtain 
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Where we have defined 



(5.10) 



Now the trajectories s and s' (s and §') arrive to the same final point r — Sf/2 (r + Sr/2). Since the initial wave-packet 
is concentrated around r , we can further simplify and work with trajectories s and s' (s and §') that have the same 
extreme points. Therefore, we have 
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By the same considerations as before, we can reduce all four trajectories to start at the center r of the initial 
wave-packet (Fig. 15) 
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FIG. 15: For fairly localized initial wave-packet, the four classical trajectories contributing to the LE can be reduced to those 
starting at its center r . 

with 

■S = P s P s < + Ps - Pg' , (5.14a) 

T = P S + Ps> -Pa-Py • (5.14b) 

Given that 
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and since the pairs of trajectories (s, s') and (s, s') have the same extreme points, the dominant contribution to /s 
will come from the terms with s — s' and s = s' . Such an identification minimizes the oscillatory phases of the 
propagators, and corresponds to the first diagonal approximation of the calculation of Sec. II and Rcf. 18. Within 
such an approximation we have 



/ E (r,i) = (j^j J ddr J2 °s C s cxp 



exp 



Sr 



(P 



2ft 2 



2^ fPs + Ps 

ft 2 



Po 



UAS s (r-°{,r ,t)-AS s (r + S {,r ,t 



(5.16) 



As in Eq. (2.10), AS s j is the extra contribution to the classical action that the trajectory s (s) acquires by effect of 
the perturbation S. 

We have two different cases, depending on whether or not there are trajectories leaving from ro with momentum 
close to po that arrive to the neighborhood of r after a time t. In the first case r is in the manifold that evolves 
classically from the initial wave-packet (Fig. 16). Such a contribution is dominated by the terms where the trajectory 
s remains close to its partner s, and calling this diagonal component, we get 
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Assuming, as in Sec. II and Ref. 18, that Ho stands for a chaotic system and that the perturbation £ represents a 
quenched disorder, upon average we obtain 
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FIG. 16: Classical trajectories in the manifold that evolves classically from ro to r, representing the diagonal component of 
/e- The action differences AS of trajectories s and s are correlated. The shaded regions depict the initial and final classical 
densities. 

where A is given by Eq. (2.30). We therefore have 
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and the corresponding contribution to the Loschmidt echo is 
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As in Eqs. (2.9) and (2.23) we have used C as the Jacobian of the transformation from r to p. Now the dominant 
trajectories are those starting from ro and momentum po. We are then back to the case of the previously discussed 
(Eqs. (2.31) and (3.36)) diagonal contribution. 



M d (t) ~ A e- xt , (5.21) 

where C = {m/t) d e~ xt is assumed, and A — [ma j A x l 2 t) d . The decay rate of the diagonal contribution is set by the 
Lyapunov exponent A, and therefore independent on the perturbation S. 

The second possibility we have to consider is the case where there does not exist any trajectory leaving from r 
with momentum close to p that arrives to the neighborhood of r after a time t. It is a property of the Wigner 
function that in the region of phase space classically inaccessible by X t the points r half-way between branches of the 
classically evolved distribution will yield the largest values of /s (Fig. 17). The trajectories s and s visit now different 
regions of the configuration space, therefore the impurity average can be done independently for each of them. As in 
Eq. (2.20), we have 
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Such an average only depends on the length L = v t of the trajectories. Thus, after average the non-diagonal term 
writes 
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FIG. 17: Non-diagonal classical contribution to the LE given by trajectories departing from ro and arriving to points equidistant 
from the point r where the Wigner function is evaluated. The action differences AS associated with both trajectories are 
uncorrelated. 



The trajectory s (s) goes between the points r and r=F Sr/2. That is why the largest values of /£ d (r,i) are attained 
when r is in the middle of two branches of the classically evolved distribution. Other points r result in much smaller 
values of /£ d (r,i), since the classical trajectories that go between r and r =p Sr/2 require initial momenta p s (p-) 
very different from p . Thus, exponentially suppressed contributions result. 
The non-diagonal contribution to the Loschmidt echo can now be written as 
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As in Eqs. (2.23) and (3.24), we have made the change of variable from r to p, and accordingly, we have obtained the 
non-diagonal contribution to the LE 18 . As discussed before, such a contribution is a Fermi Golden Rule like 22 . In the 
limit of H —* our diagonal term, Eq. (5.21), obtained from the final points who follow the classical flow, dominates 
the LE, consistently with our findings of Sec. IV B. 



C. Decoherence and emergence of classicality 



Decoherence in a quantum system arises from its interaction with an external environment, over which the observers 
have no information nor control 48 ' 65 ' 66 . The states more sensitive to decoherence are those with quantum superposi- 
tions (Schrodinger cat states), since they depend strongly on the information coded in the phase of the wave- function, 
which is blurred by the interaction with the environment. 

The studies of decoherence have traditionally considered one-dimensional systems, and often ignored the crucial role 
of its underlying classical dynamics 67 . On the other hand, it has been proposed 19 and later corroborated numerically 68 , 
that for a classically chaotic system the entropy production rate (computed from its reduced density matrix) is given 
by the Lyapunov exponent. Moreover, as shown in Ref. 18 and thoroughly discussed in this work, the decay rate 
of the Loschmidt echo in a multidimensional classically chaotic system becomes independent on the strength of the 
perturbation that breaks the time reversal between two well-defined limits (and set by the Lyapunov exponent) . The 
connection between decoherence and Loschmidt echo has been discussed in Refs. 18,69 and has induced us to note as 
Tcjy the relaxation rate of the LE. 

Decoherence is typically analyzed through the time decay of the off-diagonal matrix elements of the reduced density 
matrix (where the environmental degrees of freedom of the total density matrix of the system and its environment are 
traced out) , while the wave-function superposition defining the LE can be cast as a trace of reduced density matrices 
or Wigner functions evolving with different Hamiltonians (Eq. 1.2). Zurck has recently proposed to consider the 
relevance of sub-Planck structure (in phase-space) of the Wigner function for the study of quantum decoherence 36 . 
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Considering the example of a coherent superposition of two minimum uncertainty Gaussian wave-packets (of width 
a, centered at ±xo, and with vanishing mean momentum) in a one dimensional system, where the Wigner function 
(up to a normalization factor) is given by 
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it is clearly seen that in phase-space this distribution presents two spots located around ±xo positive with positive 
values, and between them an oscillating structure taking large positive and negative values (called interference fringes 
for their similitude with a double-slit experience). It has then been proposed that the fringes substantially enhance 
the sensitivity of the quantum state to an external perturbation. A strong coupling with an environment suppresses 
the fringes, and the resulting Wigner function becomes positive everywhere and similar to the corresponding Liouvillc 
distribution of the equivalent classical system (with statistical mixtures instead of superpositions) 66 . Jacquod and 
collaborators 41 have contested this approach, by demonstrating that the enhanced decay is described entirely by the 
classical Lyapunov exponent, and hence insensitive to the quantum interference that leads to the sub-Planck structures 
of the Wigner function. 

Working with the superposition of two Wigner functions (as in the case of the echo) and with genuinely multidi- 
mensional classically chaotic systems allows us to give a consistent description of the connection between quantum 
decoherence and the Loschmidt echo and the emergence of classical behavior. 

In the previous sections, from the semiclassical evolution of the Wigner function we were able to identify the non- 
diagonal component M nd as the contribution to the LE given by the values of the Wigner function between the 
branches of the classically evolved initial distribution (Fig. 17). Using the example of the two Gaussians of Eq. (5.25) 
(but keeping in mind that the situation is more complicated since our chaotic dynamics yields a much richer structure 
in phase-space), we see that in the region between branches both of the Wigner functions contributing to (1.2) are 
highly oscillating, and quite different from each other. The overlap, which is perfect for zero coupling (ensuring the 
unitarity requirement) is rapidly suppressed with increasing perturbation strength. As discussed earlier in the text 
(see also Refs. 18 ' 22 ), when M nd is the dominant contribution to M, we are in the Fermi Golden Rule regime. We have 
seen that this weak perturbation regime collapses as H — > (Eqs. 2.32 and 3.39). 

Beyond a critical perturbation, the diagonal component Md takes over as the dominant contribution to the LE, 
and is given by the values of the Wigner function on the regions of phase space that result from the classical evolution 
of the initial distribution. This is the Lyapunov regime, where the decay rate of M{t) is given by A. Notice that this 
behavior is still of quantum origin, as we are comparing the increase of the actions of nearby trajectories by the effect 
of a small perturbation, assuming that the classical dynamics is unchanged. The behavior in the Lyapunov regime 
does not simply follow from the classical fidelity, where the change in the classical trajectories is taken into account, 
and the finite resolution with which we follow them plays a major role. The upper value of the perturbation strength 
for observing the Lyapunov regime is a classical one, i.e. h independent (£t r — L in Sec. II D and Eq. (4.3)). 

For stronger perturbations (see discussions in Sects. II D and IV B) the classical trajectories are affected and the 
decay rate of the LE is again perturbation dependent. The Wigner function approach to the LE also helps to develop 
our intuition about the quantum to classical transition. The Lyapunov regime is the correct classical limit of a chaotic 
system weakly coupled to an external environment. 



VI. CONCLUSIONS 



In this work we have studied the decay of the Loschmidt echo in classically chaotic systems and presented evidence 
for the universality of the Lyapunov regime, where the relaxation rate becomes independent of the perturbation, 
and given by the Lyapunov exponent of the classical system. Using analytical and numerical calculations we have 
determined the range (in perturbation strength) of the Lyapunov regime, its robustness respect to the classical limit, 
the form of the perturbation, and the average conditions. 

We presented semiclassical calculations in two different Hamiltonian systems: a classically chaotic billiard perturbed 
by quenched disorder, and a Lorentz gas where the perturbation is given by an anisotropy of the mass tensor. In the 
later model, the numerical simulations were found in good agreement with the analytical calculations, and showed 
that the Lyapunov decay extends arbitrarily beyond the Ehrenfest time (where the quantum-classical correspondence 
is no longer expected to hold). 

Using a Wigner function representation, we have been able to present an alternative interpretation of the two 
contributions to the Loschmidt echo. The non-diagonal (Fermi Golden Rule) regime obtained for weak perturbation 
was shown to arise from the destruction of coherence between non-local superpositions thus destroying the non- 
classical part of the distribution. In contrast, the diagonal (Lyapunov) regime obtained for stronger perturbation 
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or more classical systems was shown to be given by the classical part of the evolved initial distribution. Thus, 
the Lyapunov regime is associated with the classical evolution (even though is of quantum origin), while the Fermi 
Golden Rule has a purely quantum nature. In this way, the persistence of the Lyapunov regime after Ehrenfest 
time is understood as the emergence of classical behavior due to the fast dcphasing of the purely quantum terms. 
This is in consistency with the understanding of the quantum-classical transition in quantum systems coupled to an 
environment driven by the decoherence 65 . 

The existence and universality of an environment-independent regime and its consequence in the phase-space be- 
havior of the Wigncr function provide a highlight on the connection between the Loschmidt echo and quantum 
decoherence. Such a connection, as well as the experiments testing the universal behavior, are promising subjects for 
future research. 

The universal behavior of the Loschmidt echo requires an underlying classically chaotic system, like the ones we 
have considered in this work. Hamiltonian systems with regular dynamics have been shown to exhibit an anomalous 
power-law for the decay of the Loschmidt echo 30 . This behavior is quite different from the one we obtain for chaotic 
systems. Therefore, we see that the Loschmidt echo constitutes a relevant concept in the study of Quantum Chaos 13 . 
Such a connection, clearly deserves further studies. 

The Loschmidt echo in the Lorentz gas has been recently calculated for short times 70 , and a rate given by twice 
the Lyapunov exponent has been proposed. It would be interesting to investigate if the difference in times scales is 
responsible for the departure from our results. 
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APPENDIX A: CLASSICAL DYNAMICS WITH AN ANISOTROPIC MASS TENSOR 

Let us assume a particle in a free space with mass tensor m surrounded by an infinite potential surface (hard wall). 
Suppose that the particle departs from a point r at time to and arrives to a final point r at time t. We must calculate 
the time t c and position r c along the surface at which the particle collides. The action along the trajectory is 



(r c - r )m(r c - r ) (r - r c )m(r - r c ) 

2(t c -to) 2(t-t c ) ' 1 1 

We can solve the problem by minimizing the action, taking the derivative of Eq. (Al) along the surface. Introducing 
unitary vector n normal to the surface at the point of collision, we can express the minimization condition as 

nxVr c S = 0. (A2) 
Denoting the initial and final velocities as Vj = (r c — r )/(t c — to) an d V/ = (r — r c )/(i — t c ), we can write 



nxm(vj-v/)=0. (A3) 

This, along with the conservation of energy E — vmv/2, results in Eqs. (3.15). The same result is obtained in the 
case of stretched boundaries. 



APPENDIX B: NUMERICAL METHOD TO SIMULATE THE QUANTUM DYNAMICS 



In order to compute the quantum dynamics of the system we resort to a lattice discretization (tight-binding model) 
in a scale a much smaller that the wavelength of the packets. This condition is relevant to accurately recover the 
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dispersion relation of the free particle Ek = h 2 k 2 /(2m), from the energy spectrum of the (open boundaries) discretized 
system, 

2h 2 h 2 

E k = k 5- (cos(fc x a) + cos(fc„a)) . (Bl) 

Setting the lattice step as the unit length (a = 1) we typically worked with R = 20, except for the calculations 
in Sec. IV B where, in order to keep the precision for smaller wavelengths, a was reduced keeping the product ka 
constant. 

The discretization results in a Hamiltonian matrix whose diagonal elements are the on-site energies. The off-diagonal 
elements are hopping terms V = ft 2 / (2ma 2 ) which is the maximum kinetic energy represented by the discretization. 

The quantum dynamics on the lattice was carried out using a Trotter-Suzuki algorithm 71 , which is a remarkably 
precise and efficient numerical method. At the lowest order, it is a decomposition of the evolution operator U for a 
small time r in a product of analytically solvable evolution operators. Typically one searches for a way to write the 
Hamiltonian of the system as H = J2k ^fe> where Ttk are 2x2 matrices, and thus 

Q 

U(t) = cxp pir/h] ~ U(t) = Y[ exp [iW fc r/fi] , (B2) 

k 

where Uk(r) = exp[iHfc-r/fi.] are rotation matrices. 

The highest orders involve a fractal decomposition of t that preserves the unitarity of the approximated evolution 
operator. In our calculations, a fourth-order algorithm with a time step r = O.lh/V was precise enough for the time 
regime of interest. 
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